Fully coupled forward-backward stochastic 
dynamics and functional differential systems 

Matteo Casserini* and Gechun Liang^ 



Department of Mathematics, ETH Zurich, Switzerland, 
matteo . casseriniSmath . ethz . ch 

^Oxford-Man Institute, University of Oxford, U.K., 
gechun . liangOoxf ord-man . ox . ac . uk 



December 22, 2011 



Abstract 

This article introduces and solves a general class of fully coupled 
forward-backward stochastic dynamics by investigating the associated sys- 
tem of functional differential equations. As a consequence, we are able 
to solve many different types of forward-backward stochastic differential 
equations (FBSDEs) that do not fit in the classical setting. In our ap- 
proach, the equations are running in the same time direction rather than 
in a forward and backward way, and the confiicting nature of the structure 
of FBSDEs is therefore avoided. 
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1 Introduction 



Due to their central role at the intersection between stochastic analysis, mathe- 
matical finance and partial differential equations, backward stochastic differen- 
tial equations (BSDEs) and forward-backward stochastic differential equations 
(FBSDEs) have been subject of extensive research during the last two decades. 
While linear BSDEs had already been introduced by Bismut [2] in 1973, it was 
only after the seminal work of Pardoux and Peng [19] in 1990, who first studied 
the general non-linear case, that BSDEs gained considerable attention. Since 
then, the importance of the theory of BSDEs increased dramatically, finding 
numerous applications in stochastic control theory, PDE theory, mathematical 
finance and many other fields. We refer the reader to the books [10, 18] and the 
surveys [9, 11] for an extensive overview of BSDEs and their applications. 

While simple types of decoupled FBSDEs had already been considered by 
Pardoux and Peng [20], the study of fully coupled FBSDEs has been initiated 
by Antonelli [1], who proved the existence and uniqueness of local solutions. 
The solvability of fully coupled FBSDEs was later studied by several authors, 
and mainly three types of methods have been proposed so far, each having 
its constraints and which do not cover each other. The first is the method of 
contraction mapping, introduced in the local case by Antonelli [1]: his work 
was later developed by Pardoux and Tang [21] to solve FBSDEs globally under 
additional monotonicity conditions. Later on, motivated by the method of con- 
tinuation in PDE theory, Hu and Peng [13], Peng and Wu [22] and Yong [24] 
solved FBSDEs on arbitrary intervals by relying on a different type of mono- 
tonicity assumptions on the coefficients. Finally, Ma et al. [16] introduced the 
well-known four-step scheme, which links FBSDEs and quasilinear PDEs: in 
this case, the coefficients have to be deterministic and satisfy strong regularity 
assumptions. This method was further developed by Dclaruc [4], who relaxed 
the regularity assiimption on the coefficients by combining the four-step scheme 
with the contraction method. In the last years, Zhang [25,26] and more recently 
Ma et al. [17] developed a so-called decoupling scheme to solve FBSDEs with 
random coefficients, which unifies most of the existing results at least in the 
one-dimensional case. For an extensive account of FBSDEs, we refer to the 
book by Ma and Yong [18]. 

The purpose of this article is to introduce a general class of fully coupled 
forward-backward stochastic dynamics on a general filtered probability space, 
which contains classical FBSDEs as a special case. Inspired by the recent work 
on Lipschitz BSDEs of Liang et al. [15] , we study the solvability of these forward- 
backward dynamics by introducing an appropriate functional differential system 
of the form 

dXt = n{t, Xt, y{X, V)t, (Mix, V))t)dt + a{t, Xt,yiX, V)t,C\M{X, V))t)dWt, 

dVt = fit, Xt,y{x, v)t,/:HM{x, v))t)dt, 

Xo =x, Vo= 0, 

where fi, a, f are random functions, £^ are general abstract operators, 

and M., y are given, for a random function cj), by 

M{X, V)t = M't'iX, V)t := E[cj>{XT) + VrlJ't], 
y{x, v)t = y (X, V)t := M'^iX, V)t - Vt. 
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This method is purely probabihstic and does not rely a priori on the existence 
of martingale representations. 

In particular, our approach shows that FBSDEs can be reformulated as 
functional differential equations defined in a forward way: since both equations 
are running in the same time direction, this avoids the conflicting nature of the 
structure of FBSDEs. More important, our results allow us to consider a more 
general class of forward-backward stochastic dynamics which are beyond the 
existing framework, and extend the results of Liang et al. [15] from backward 
systems to fully coupled forward-backward systems. 

After introducing an appropriate framework and defining properly the prob- 
lem, we study its local solvability and derive our main result: the existence of a 
unique local solution to the functional differential system under Lipschitz and 
monotonicity conditions on the coefficients and under specific Lipschitz assump- 
tions on the operators £*. In particular, the conditions on £* are rather mild and 
allow to consider many types of operators different from the usual martingale 
integrand processes treated in classical FBSDEs: as a consequence, we can solve 
within our framework many different types of forward-backward equations that 
do not fit in the classical FBSDE setting. To emphasize the generality of these 
assumptions, we present several examples of possible operators and potential 
financial applications. 

In the second part of the article, we discuss the solvability of the system on 
arbitrarily large time intervals. This is however more problematic: indeed, it 
appears impossible to study such an extension without defining the operators O 
explicitly, and one has to consider the problem separately for each choice of 
We conclude the article by presenting a study of the case where the filtration is 
Brownian and the operators £* are given by Ito's representation. 

The article is organized as follows. First of all, to provide some intuition, 
we give in Section 2 a brief overview of the functional differential approach in a 
simple Brownian setting. In Section 3, after introducing a more general frame- 
work, we give a rigorous definition of our class of forward-backward dynamics 
and the associated functional differential system. Section 4 is then dedicated 
to the existence and uniqueness of solutions to the latter system for sufficiently 
small time horizons. Finally, in Section 5 we discuss the general problem of ex- 
tending the solution to arbitrarily large time intervals, and study in particular 
the case of classical Brownian FBSDEs. 

2 The functional differential approach 

We would like to begin by providing the reader with some intuition of the 
approach we are going to use in the sequel, which is inspired by the work of 
Liang et al. [15]. We first present the following elementary, but very illustrative 
result derived in [15]: 

RerriEtrk 2.1 (Liang et al. [15]). For T > 0, assume that we have a special 
semimartingale {Xt)t<^[o,T] on a probability space (fi, J", {Tt)t^[o.T]i P) satisfying 
the usual assumptions, and let the terminal value Ft = C G L^{^t) be given. 
Furthermore, assume that the canonical decomposition of Y is given by 

Yt=Mt- Vt, 
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where M is a martingale and V a predictable process of finite variation with 
Vo = 0. Then, if Vr is integrable, it is easy to verify that, for all t G [0, T], 

Mt = E[MT\Tt] - + VrlTt], 
Yt = Mt-Vt^ + VT\Ft] - Vt. 

In other words, the semimartingale Y and the martingale M can be expressed 
as operators of the terminal value ^ and the finite variation process V. 

We show now, with the help of some intuitive arguments, how this remark 
can lead us to an alternative formulation of the classical FBSDE problem. Let 
(n, J^, P) be for the moment a complete probability space with an m-dimensional 
Brownian motion W = {Wt)te[o,T] and the corresponding augmented filtration 
{J-'t)te[o,T], and consider a classical fully coupled FBSDE of the form 

'dXt = Xt,Yt, Zt)dt + a{t, Xt, Yt, Zt)dWt, 
< dYt = -fit, Xt, Yt, Zt)dt + ZtdWt, (2.1) 

Xo = X, YT = (t>{XT), 

where the functions ^ : x [0, T] x M" x R** x R'*^™ ^ R", ct : x [0, T] x K" x 
X Kdxm jgnxm^ / : X [0, T] X M" X R"^ X R''^''" -i^ R'', : n X R" -)> R"^ 

satisfy the usual measurability and integr ability conditions. 

Assmue now that the above FBSDE has a solution {X, Y, Z). Since is 

the terminal value of the semimartingale Y, Remark 2.1 induces us to introduce 

the operators M'^ and by defining 

M'^iX, V)t := E[cj>{XT) + VT\Tt], 
y^{X,V)t:=M'^{X,V)t-Vt, t€[0,T], 

for any processes X, V such that G L^{Ft), Vt G L^{J-t)- Then, by 

Remark 2.1 and the definition of the operators and M'^ , it seems plausible 
to associate the above FBSDE to the following system of forward functional 
differential equations 

'dXt = n{t, Xt, y^ix, v)t,z^{x, v)t)dt + a{t, Xt, y^{x, v)t,z^{x, v)t)dWt, 
< dVt = f{t, Xt,yHx, v)t, z^ix, v)t)dt, 

Xo = x, Vo = 0, 
where Z'^ is given implicitly via Ito's representation theorem by 

M'f'iX, V)t = E[M'>'{X, V)t] + I Z't'iX, y),rfW«. 

Jo 

The peculiarity of these stochastic differential equations consists in the fact 
that the coefficients /i, cr, / depend not only on the behaviour of the solution 
process {X, V) up to the present value, but also on the terminal value {Xt, Vt) 
of the solution: such stochastic differential equations are not standard, and 
for this reason they arc called functional differential equations (note that the 
term "functional differential equations" is often used in the literature to refer 
to stochastic delay differential equations: however, contrary to the latter, the 
drivers of the equations studied here do not have any delay in the past, but 
rather in the future). 
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3 Fully coupled functional differential systems 



In this section, we will show how the approach presented above can be made 
rigorous and extended to a much more general framework. To this end, we first 
need to introduce some notation: in the following, we fix T > 0, and assume 
that we are given a complete probability space (O, P) together with a general 
filtration F = {J^t)telo.T] satisfying the usual assumptions. We recall that every 
F-martingale has under these conditions a cadlag version, which we will always 
choose. Moreover, we denote by V the predictable a-field with respect to F and 
assume that an m-dimensional Brownian motion W = {Wt)t£[o t] is defined on 
(r!,J-,F,P). 

For fc, / e N, I • I denotes the Euclidean norm on M'^, respectively the Hilbert- 
Schmidt norm on R*^^', and R'^^' will often be identified with R'^'K We define 
^2([o,T],M'') as the space of all processes V : n x [0,T] ^ R"^ continuous and 
adapted such that = and £[suptg[o_T] < oo, while ^^([0,T],R'^) 

denotes the space of all square integrable R'^-valucd martingales on [0,T]. Both 
y^{[0,T],R'^) and ^^{[0,T],R'^) are endowed with the norm 



sup \Vt\ 
te[o,T] 



and note that (^^([0, T], R''), || • |1,s^2[o_t]) is then a Banach space. Sometimes, 
we will also need the direct sum space ( [0, T] , R'') © qq, T] , M'') , endowed 
with the same norm || • ||^2[o,t]- 

In the following, we will consider a particular generalization of the FBSDE 
(2.1) on the filtered probability space (fi, J^, F, P). More exactly, we will assume 
that fjL, u and / depend on some general processes jC'^{M), i = 1,2,3, instead 
than on Z, where is for each i an abstract operator defined on ./#^([0, T], K*^) 
and taking values, for some pi G N, in the space of p^-dimensional adapted 
processes (the codomains of will be further specified later). As we will see 
later more in detail, this substitution will allow both to take into account the 
generality of the filtration F and to treat locally other types of forward-backward 
equations not fitting in the classical framework. This generalization of the 
FBSDE (2.1) leads us to the following system: 

dXt = n{t, Xt,Yt, C^{M)t)dt + a{t, Xt, Yt, C\M)t)dWt, 
dVt = -fit,Xt,Yt,CHM)t)dt + dMt, (3.1) 
Xo = x, Ft = </>(Xt), 

where /x : x [0, T] x R" x R'' x R^i ^ R", cr : x [0, T] x R" x R'' x Rp^ 
R"^™, f : n X [0,T] X R" xR'^ X RP^ ^ R'', (f) : n X R" ->■ R'* have the 
necessary measurability and integrability properties. Note that the system is 
then completely determined by (/i, a, f, (j), L"^ ,£? ,£?). 

Definition 3.1. A solution to (3.1) is a triplet of processes {X,Y,M) such that 
X e y^{[0,T],R"), Y e y^([0,T],R'^) ® ^^([0,T],R'^), M G ^^{[0,T],R'^), 
and satisfying the integral formulation of (3.1). 

We will call such systems fully coupled forward-backward stochastic dynamics. 
As mentioned in the previous section, a viable approach to study the solvabil- 
ity of this system is to reformulate the problem with the help of appropriate 
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functional differential equations. To this end, wo denote by "^^^ the class of 
M"-valued adapted processes X on [0, T] such that (t>{XT) G L^{Tt), by the 
class of M'^-valued adapted processes V on [0,T] such that Vt S L^[Ft), and 
we define the operators M'^ and on x by 

M-^IX, V)t := i^[</.(XT) + Vrl^t], 

y (X, y)t := A^'^(x, y)* -Vu te [o, t]. 

For the rest of this article, we will drop the dependence of Al'* and on 
by writing Ai and y. With the help of the operators Ai and 3^, we reformulate 
the problem (3.1) as the following fully coupled system of forward functional 

differential equations: 

( dXt = nit, Xt,y{X, V)t,C' {Mix, V))t)dt + ait, Xt, 3^(X, V)t, £:'{M{X, V))t)dWt, 
I dVt = fit, Xt,yiX, V)t,jC^iMiX, V))t)dt, (3.3) 
[xo=x, Vo = 0. 

Definition 3.2. A solution to (3.3) is a pair of processes (X, V) such that 
iX,V) G ^2Qo,T],M") X S^^i[0,T],W^) and satisfying the integral formulation 
of (3.3). 

Such systems will be called fully coupled functional differential systems. It 
is not diSicult to show the equivalence of the systems (3.1) and (3.3): 

Lemma 3.3. The fully coupled forward-backward system (3.1) has a solution if 
and only if the functional differential system (3.3) does. 

Proof. If iX, Y, M) solves (3.1), then we obtain a solution of (3.3) via the canon- 
ical decomposition of the semimartingale Y. Conversely, if (X, V) is a solution of 
the functional diff'erential system, then {X, y{X, V),MiX, V)) solves (3.1). □ 

An immediate observation is that the functional differential system (3.3) 
has a more homogeneous structure than that of the original problem. Indeed, 
while the forward-backward dynamics (3.1) consist of a forward and a backward 
equation of different nature, both the functional differential equations in (3.3) 
are running forward in time and show a similar dependence of the coefficients on 
both the present and the terminal values of the solution processes. In particular, 
this homogeneity allows to rewrite the problem more compactly as 

'dUt = ^{t,7riiUt),yiU)t,CHMiU))c,CHMiU)),)dt 

+Eit,MUt),yiu)t,c^iMiu))t)dWt, 

Uo = ix,0)^, 

where U = iX, V)"^ , tti : — > R" is the projection on the first n compo- 
nents, S = (cr,0)'^, and for t € [0,T], (.x, y, zi, zg) G K" x ]R<^ x Rp^ x Rp^, 
^'(f, .T, y, zi. z;j) = ifiit,x,y, zi), fit,x,y, z^))"^ . For the rest of the article, we 
prefer however to consider the system as formulated in (3.3), since it will be 
more convenient to treat the coupling between X and V in such a framework. 
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4 Existence and uniqueness of local solutions 

Because of Lemma 3.3, we will now focus our attention on the existence and 

uniqueness of solutions to the fully coupled functional differential system (3.3). 
As a first step, we study the local solvability of the problem and introduce 
sufficient monotonicity and Lipschitz assumptions on the coefficients n, a, f, (j) 
and the operators £^ and C'^ . 

In the following, we denote by ^^([0,T],M') the space of P-measurable 

processes : x [0,T] such that ||iI||^2[o,T] ■= E[Jq \Ht\^dt] < oo, 

where I € N. First of all, we shall assume that the coefficients fi, a, f, (p satisfy 
the following assumption: 

Assumption (Al): The functions ^ : x [0, T] x M" x R'* x M", 
a ■.nx[0,T]xR" xR"^ X Rp^ ->■ R"><™, / : x [0, T] x M" x M'* x M^^ _^ 
and : f2 x K" ^ R** satisfy Assumption (Al) if there exists a constant C > 
such that: 

(Al.l) For any {x,y, zi, Z2, z^) e R" x R'' x Rp^ x Rp^ x Rp^ , the processes 
^{■,x,y, zi), a{-,x,y,Z2) and f{-,x,y,zz) are V -measurable and (p{x) 
is -measurable. 

(A1.2) For every {x,y, zi), {x' ,y' , z[) G R" x R'^ x Rp\ 

{x - x')^{fj.{-,x,y,zi) - n{-,x',y,zi)) < C\x - x'p, 
\l^{-,x,y,zi) - ii{-,x,y',z[)\ < C\\y - y'\ + \zi - z[\), 

X, 0, 0)1 < C(l + |a;|) dP (g) dt-a.s., 

and the function x x, y, zi) is dP ® dt-a.s. continuous. 

(A1.3) /(•, 0,0,0) € ^^{[Q,T],R'^), c7(-, 0,0,0) G jr2([0,r],R"><'") and 
4>{0) e L2(n,Md). 

(A1.4) For every (x, y, Z2, Z3), {x', y', z'^, 4) G R" x R"^ x Rp= x Rp^, 

\a{;X,y,Z2) - a(.,x',y',4)|2 < C{\x - + \y - yf + \z2 - 41'), 
\f{;X,y,Z3) - fi;x',y',z',)\ < C{\x-x'\ + \y-y'\ + \zs - 4|), 
\^{x)-cp{x')\<C\x-x'\ dP® dt-a.s.. 

Remark 4.1. The above conditions on /i, cr, / and (/) arc quite standard in the 
theory of FBSDEs (see for instance [18]). Moreover, the reader can easily verify 
that the condition (A1.2) could be replaced by the following stronger, but more 
standard assumption: 

(A1.2') /x(-,0,0,0)eJf2([0,T],R") and, for {x,y,zi),{x',y',z[)&R''xR'^xRP\ 
\fi{t, x, y, zi)-fj,{t, x', y', z[)\ <C{\x-x'\-\-\y-y'\-\-\zi-z[\) dP^dt-a.s.. 

In particular, the assumptions on (j) allow us to derive the following Lipschitz 
estimates on the operators y and M., which will be essential in the sequel: 
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Lemma 4.2. Assume that (f> satisfies the conditions in (Al). Then we have 
that, for the operators introduced in (3.2), 

Moreover, for any X,X' e =y2([o,T],M") and V,V' € ^^{[0,T],m^), 

\\y{x,v) - y{x',v')\\^2[o,T] < 2C\\x - x'\\j..^o^T]+nv -v'\\^2[o,t], 

\\M{X, V) - M{X', y')IU=[o,T] < 2C\\X - X'||^2[o,T] + 2||y - V'\\ 
Proof. We first prove the second assertion. By tlie triangle inequality, 

\\M{X, V)-M{X',V')\\.y2[„^T] 

< \\EWXt) - <^(X^)|J-.]||^,[,^^] + \\E[Vt - vm\\^2io,T] 
= E 



f 


\E[Vt - 






1/2 




|2' 


+ E 


sup 









sup \E[^{Xt) - (^{X'^)\Tt\^ ' " 



-te[o,T] 

and therefore, by Doob's inequality and the assumption on cj), 
||Al(X,F)-Al(X',yOII^^[o,r]<2(i?[|0(XT)-0raP]'/+i?[|VT-l4P]'^') 

<2C\\X - X'||^2[o,T] + m - y'||^2[o,T]- 

The estimate for |13^(X, F) — 3^(^', V)||3^2[o,t] then follows by the application 

of the triangle inequality. □ 

The next step consists in introducing appropriate conditions for the abstract 
operators £^ and C^. These are given by the following Lipschitz and bound- 
edness assumptions, which are the same as introduced by Liang et al. [15]. 

Assumption (LI): The operator C satisfies Assumption (LI) if: 

(Ll.l) C maps .y^^{[Q,T],W^) into 0^{[Q,T],W), where ^^([o, T], Rp) is ei- 
ther the space M''^{[0,T],W) or y^{[Q,T],W). 

(LI. 2) C is hounded and Lipschitz continuous, i.e. there exists a constant 
K>Q independent of T such that, for all M, M' e ^^([0, T], I 

ll^(M)|U2[o,T] <if||M||^2[o,T], 

||£(M) - £(M')||^2[o,T] < K\\M - M'||^2[o,T]. 



We will see in Theorem 4.5 that Assumption (LI) is enough to guarantee 
the local solvability of our system without a concrete specification of the op- 
erators In particular, as anticipated in the Introduction, the weakness of 
Assumption (LI) allows to consider many different types of operators within our 
framework: we emphasize its generality by giving several examples of possible 
operators, and we present potential financial applications. 

Examples 4.3. We begin with the classical case of integrand processes gener- 
ated by martingale representations. 
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(i) Assume that (^t)te[o,T] is the augmented filtration generated by the Brow- 
nian motion W, and take ff^i[0,T],W) = ^^^([0, T], R'^^™). Then, we 
define jC : ^'^i[0,T],R'^) -)> ^^^([0, T], implicitly via Ito's repre- 
sentation theorem by 

Mt = Mo+ I L{M),dWs, t G [0, T]. 
Jo 



By Ito's isometry we have that 



\\c{m)\\1,.,,,,,.=e' 



\C{M)t\''dt 







= E[{Mt - Mo)^] = E[Ml] - E[Ml] < ||M||,2^.[o,t], 

and the Lipschitz property follows by the linearity of £. Note that in the 
case where JC^ = C? = C? = C, the system (3.1) is reduced to a classical 
FBSDE: this shows that classical FBSDEs can be seen as a special case 
of the forward-backward stochastic dynamics (3.1). 

(n) Let now (.^t)te[o,T] be a general filtration with just the usual assumptions. 
As in the previous case, we can take e'^{[{),T],W) = Jf ^([0, T], M''^'"), 
and define C : .^^{[O.T].R'^) .9if'^{[0,T],R'^''"') implicitly via the or- 
thogonal decomposition with respect to W, i.e. 

/■* 



= / C{M)sdWs+Nt, t€[0,T], 
Jo 



where A'' is some martingale orthogonal with respect to W. The reader 
may easily notice the connection between this operator and generalized 
BSDEs (see [10]), and we can prove similarly to Example 4.3 (i) that jC 
satisfies (LI), by applying the Burkholder-Davis-Gundy inequality instead 
of Ito's isometry. 

We can thus study generalized fully coupled FBSDEs within our frame- 
work, and the choice of the filtration allows us to consider, in typical 
financial applications, the case of incomplete markets. An illustrative ex- 
ample is that of a large investor trading in an incomplete market: since this 
investor buys and sells large amounts of assets, it is reasonable to assume 
that his trading strategy affects the prices of the stocks. By considering 
the corresponding hedging problem, we thus obtain a fully coupled system, 
and the incompleteness of the market leads to a generalized fully coupled 
FBSDE. For more details, we refer the reader to [18] . 

Examples 4.4. While martingale integrand processes are the case most studied 

in the literature, they are not the only class of operators fitting in our framework: 
there are indeed several other classes of non-local operators, not considered in 
the classical FBSDE literature, which satisfy Assumption (LI). Let us give 
some examples. 

(z) Assume that {J^t)telo,T] just satisfies the usual assumptions. We take 
^2([0,T],RP) = ^2([o,T],R'^), and £ : ^^{[0,T],R'^) S^^{[0,T],R'^) 
is simply defined by £(M) := M; in this case. Assumption (LI) becomes 
trivial. In a financial context, £(M) may represent the risky part of the 
claim Y = M -V. 
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(n) Let for simplicity d = 1. Fix T > and assume that, for T < T, 
F = (^i)te[o,T] is such that all martingales with respect to F are contin- 
uous. Let ^2([o,r],RP) = jr2([0,T],M), and define £ : ^^i[0,T],R) 
^2([o,r],M)by 



£(M)i := ^E[{M)t^T\J't], M G^^i[0,T],R), te [0,T], 

where for notational simplicity {M)t^T '■= {M)t — {M)t- Then, by the 
Kunita-Watanabe and the conditional Cauchy-Schwarz inequalities, we 
have that 



E [{M, M\t I J^t] < E [I (M, M')t,T\ \J't]<E y {M)t,T{M')t,T 



< jE[{M)t,T\J'tWE[{M')t,T\J't], 



and therefore, by the bilinearity of {■)t,T, 

|£(M)t - £(M')t|' = E[{M)t,T\J't] + EliM^TlJ't] 



- 2^ E[{M)t,T\Tt]yJ E[{M')t,T\H 
< E[{M - M\T\Tt] 



for all t > 0. Hence, by Fubini's theorem 

||£(M)-£(M')|l^.[o,Ti=£^' 



<E 
= E 



C{M)t - C{M')t\ dt 

I E[{M-M\T\Tt]dt 
Jo 

[ {M-M')t,Tdt 
Jo 



<TE[{M-M')t]. 



By applying the Burkholder-Davis-Gundy inequality, we finally get the 
desired Lipschitz property. The boundedness condition is obtained via 

similar computations. 

(m) We choose again for simphcity d= 1, and we assume that F = {Tt)te[o,T] 
is as in Example 4.4 (ii). We can define £ by taking &'^{[Q,T],W) = 
^2([o,T],M), and 



£ : .^2([0,T],K) ^ ^2([0,T],M), £(M)t := ^/{mYu t G [0,T]. 
By the Burkholder-Davis-Gundy inequality, we have that 



||£(M)|| 



E 



sup (M)t = E[{M)t] < K\\M\\^^,,o^^y 



On the other hand, by applying the Kunita-Watanabe inequality. 



^/{M)t- V{M')t = {M)t + {M')i-2^{M)t{M')t 

< {M)t + {M')t-2\{M,M')t\ 

< {M)t + {M')t - 2(A/, A/')i = {M - M')t, 
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and the Lipschitz property then follows by applying the Burkholdcr-Davis- 
Gundy inequality as above. We restrict for a moment to the Brownian set- 
ting to give a financial interpretation: in this case, £(M) can be explicitly 

rewritten as C{M)t = \Zs\'^ds, where Z is the martingale integrand 
in the Ito representation of M. In the usual BSDE framework for hedging 
(see for instance [11]), £(M) is then closely connected to the accumulated 
cost of the portfolio strategy: this could allow us, for instance, to consider 
storage problems within our setting. 

(iv) Wo modify the previous example by combining it with orthogonal decom- 
positions. Let d — I and F = {J't)te[o.T] &s in Example 4.4 (ii). Fix 
a martingale M, and define 7^ : ^^{[0,T],R) ^^{[0,T],R) as the 
orthogonal term in the orthogonal decomposition with respect to M, i.e. 

Mt = {Mt-n{M)t)+n{M)u t€[0,T], 

where TZ{M) is orthogonal with respect to M and Mt-U{M)t = J* Z^dMs 
for some process Z. Then, we define the operator C by 



£ : ^2([o,T],M) ^ =y2([0,T],M), C{M)t := ^{n{M))t, t € [0,T]. 
Because of the orthogonality, it is easy to check that, for all t>0, 

{M)t = {M -n{M) + n{M))t 

= {M- n{M))t + {n{M))t > {n{M))t, 

and similarly for M—M'. Therefore, £ satisfies Assumption (LI) because 
of Example 4.4 [iii). This operator may have interesting applications in 
mathematical finance: namely, in the typical BSDE framework for hedging 
in incomplete markets, the operator TZ{M) represents the non-hedgeable 
part of the claim. While we cannot hedge this risk, we can partly incor- 
porate its effect on the price process since the coefficients fi, u and / are 
allowed to depend on ^J {TZ{M))f 

{v) As a final example, we introduce an operator intimately related to back- 
ward equations with time delayed generators (see [7]). Let T > be 
fixed. For T < T, let {J-'t)te[o,T] satisfy the usual assumptions, and 

^2([0,T],RP) = ^2(-[o,T],R'^^'"). Motivated by the framework intro- 
duced by Dos Reis et al. [8], we can then define £ : ^^([0, T], R'') 
Jf'2([o,T],M'^x™) by 

C{M)t := j"" C{M)t+,az{ds), 

where £ is the operator introduced in Example 4.3 [ii), and az is a 
non-random finite measure with support in [— T, 0]. We observe that £ 
is closely related to the operator of Example 4.4 [iii) when az is the 
Lebesgue measure restricted to [— T, 0]. Moreover, by applying the change 
of integration order proved in [8] , we can show that 

1/2 



||£(M)||^2[o,T] = E 



C{M)t+saz{ds) 



dt 



(4.1) 



< az([-r,0])||£(M)|U2[o,T] < if||M||^2[o,T] 
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for some constant K — K{T), where the last inequahty foUows from the 
boundedness of £. Assumption (LI) then follows by linearity: this will 
therefore allow us to consider a special class of fully coupled forward- 
backward equations with delayed generators. 

However, the treatment of backward stochastic delayed equations in their 
full generality would require to replace y{X, V) in (3.3) by a new opera- 
tor Y(X,y), := J\y{X,V),+saYids), where ay is a non-random finite 
measure: as pointed out in [8], this kind of equation is very difficult to 
study even in the simple decoupled case, and it will not be treated here. 
For an overview of the several applications of backward equations with 
time delay to mathematical finance, the reader may consult for instance 
the article by Belong [6]. 

The broad generality of the above examples should help us understand the 
importance of the next theorem. This is the main result of this section, and gives 
the existence of a unique square-integrable solution to the functional differential 
system (3.3) on sufficiently small intervals [0, T], provided that the operators 
satisfy Assumption (LI). This gives us a great flexibility in the choice of 
allowing to study many different types of fully coupled, non-classical forward- 
backward stochastic dynamics. 

Theorem 4.5. Let fi, a, f and (p satisfy Assvmption (Al) luith respect to 
the constant C . Furthermore, assume that , £ ' satisfy Assumption (LI) and 
thatC^ satisfies (LI) with respect to ^^([o, T], Rp^) = .y^{[0,T],W^), denoting 
by K the common Lipschitz constant of , , C? . Then there is a constant 
I = £{C, K) depending only on C and K so that, for T < I, (3.3) admits a 
unique solution {X,V) in y^{[Q,T],W) x y'^{[0,T],R'^). 

Proof. Prom now on, we write || • ||^ for the norm || • ||^2[o^t] for notational sim- 
plicity. Moreover, we denote the product space J^^([0, T], ffi") x y^{[0,T],R'^) 
by =5^^ X and endow it with the norm 

\\{X,V)\\^2^^. := ^J\\X\\%, + \\V\\%,, {X,V) e ^1 X y^. 

The mapping L : x ^ x y^, 1L{X,V) := {X,V), is defined as 
follows: first, X is constructed as the unique solution in y'^{[0,T],W^) to the 
forward stochastic differential equation 

dXt = fj.{t, Xt,y{X, V)t, C\M{X, V))t)dt+a{t, Xt,y{X, V)t, C\M{X, V))t)dWt, 
Xo = X. (4.2) 

Then, once X has been obtained, V is given explicitly by the expression 

Vt= f fis,X,,y{X,V)s,CHM{X,V))s)ds, (4.3) 
Jo 

We show that the mapping L is well defined and maps y^ x yy into itself. 
First of all, we note that the existence of a unique solution in J^^([0,T],M") 
to (4.2) follows by Assumption (Al) and the results on stochastic differential 
equations with monotonous coefficients obtained, for instance, by Rozovsky [23]: 
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indeed, by setting fl{t,x) := ^{t,x,y{X,V)t, C^{M{X,V))t) and d{t,x) := 
a{t,x,y{X,V)t,C'^{A4{X,V))t), the reader can easily verify that all the con- 
ditions of [23] are satisfied, as ||3^(X, y)||^2 and ||A^(X, F)||^2 are finite by 
Lenuna 4.2. On the other hand, it is easy to verify that V e ^'^{[0,T],M.^), as 
for /o :=/(•, 0,0,0), 



^2 



< 



Vf(^\\fo\U2 + \\f{.,x.,y{x,v).,cHM{x,v)).) - fo\U^) 



<Vf(\\fo\\.^2 + aVT W l){\\X\\,r2 + \\y{X,V)l^2 + \\C^{M{X, V))\\ff2)^ 
<VT(||/o|U2+C(VTVl)(||X|U2 + p(X,y)|U2+i^||M(X,F)|U2))<oo. 

Since the pair (X, V) is a solution of (3.3) if and only if it is a fixed point 
of L, it suffices to prove that L is a contraction on x S^y for small enough 
T > 0. Let L(X\ V^) = {X\V^), L{X^, V^) = {X^, V"^), and assume without 
loss of generality that T < 1. By Ito's formula, we can compute that 

d\Xl - Xff = 2{Xl - Xffd{X' - X% + d{X' - X\ 
= 2{Xl-X^,f{^^{t,Xly{X\V')t,C\M{X\V^))t) 

- xl y{x\ v%,cHm(x^ v^))t))dt 

+ 2{Xl - X^f {a{t,Xly{X\V%X\M{X\v'))t) 

-a{t,Xly{X\V\,£.\M{X\V^))t)dWt) 
+ \a{t, xl,y{x\ V')t, C\M{X\V^))t) 

- a{t, xl y{X\ V^)u jC,\M{X\ V^))t) I'rfi 

for any t > Q. Thus, by applying Assumption (A1.2), (A1.4) and classical 
inequalities we obtain that, for a constant 6\ = 6\{C) depending only on C, 



\\X^-XY^.=E 



sup 

te[o,T] 



< 6*1 \E 
+ E 
+ E 
+ E 
+ E 



+ E 



'\xl-x^\\yix\v%-y{x^vWds 



f\Xl-Xl\'ds 
Jo 

r \Xl - X^\\jC\M{X\V^))s - /:\M{X\ V^))s\ds 
.Jo 

[ \y{x\v'),-y{x',vW^ds 

Jo 

\C^MiX\ V')), - C^iMiX\ V^)).s\^ds 



\xl - x^fiixl - x^f+\y{x\ V'), - y{x\ 

+ \C\M{X\v'))s - CHM{X', V')U^)ds' 



LJO 



(4.4) 



1/2. 
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The next step consists in deriving estimates for all the terms on the right 
hand side of this inequality. The first term can simply be estimated by 



E 



Jo 



< T\\X^-X^\\%2. 



(4.5) 



For the second term, we obtain by the Cauchy-Schwarz inequality that 



E 



Jo 



xi\\yix\v')s-y{x\v'),\ds 



< 



Because of assumption (A1.4) on cj), we can apply Lemma 4.2 and get that 



E 



\xl -x^\\yix\v^)s-y{x\v^)Ms 

< 0,t{\\x' - xY^. + \\v' - vY^^ + \\x' - xY^^) , (4.6) 

for some constant 62 = O2 (C) depending on C. To estimate the third term, one 
can apply again the Cauchy-Schwarz inequality to get that 



E 



Jo 



Xl\\C\M{X\V^))s - C\M{X\ V^))s\ds 



T 
2e' 



<^\\X^-XY^^ + l,E 



Jo 



< ^{\\X' - XY^^ + \\C\MiX\ V')) - C\MiX', v')Y 



where the last inequality is obtained by taking e = \fT if 0"^ = and e = 1 
if G'^ = S^"^ (remember that T < 1). On the other hand, by Assumption (LI) 
we know that 

\\C}{M(X\V^)) ~ C\M{X\ V^)Ye^ < K^M{X\ V') - M{X^ V^Y^., 
and by applying Lemma 4.2, 

e\ r \Xl - X^\\C'{M{X\V')), - C\M{X',V^)Uds 
.Jo 

< e3VT(\\x^ - xY^2 + \\v^ - vY^2 + \\x^ - xY^^), (4.7) 

for a constant ^3 = 0z{C,K) depending only on C and K. The fourth term is 
estimated by using the same argument as for (4.6), obtaining that 



E 



\y{xW),-y{x\v%\'ds 



< 



r2\\2 



\X'-X 



2||2 



(4.6 
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for some constant 64 = 94(C). For the fifth term, we have that 



E 



Jo 



< K^T\\M{X\V^) - M{X^,V^)\\%2, 



and hence, by Lemma 4.2, 

e' 



[ \C\M{X\V')), - C^M{X', V^)U^ds 
Jo 

< ^5^(11^^ - ^^11^2 + \\X' - XY^^), (4.9) 

for 6*5 = 6*5 (iC). It only remains to estimate the last term: for notational 
simplicity, we introduce the process A' := (y{X',v^),c^{M{x\v^))^ for i = 
1,2. By applying the Cauchy-Schwarz inequality, we can verify that 



E 



< E 



T ^ ^ ^ ^ \ 1/2- 

\Xl-Xl\Wxl-Xl\'' + \Al-Al\^)d8 



( sup \Xl-XlA"'"([ {\xl-Xl\- + \Al-Al\^)ds 



1/2- 



Vt 



sup \Xl-Xf\' 

t£[0,T] 



2Vt 



E 



and it is not difiicult to check that 



E 



[ (\Xl-X^\' + \Al-Al\')ds 
Jo 



< T 



which finally leads us to 



E 



{\Xl-X^f + \Al-Alf)ds 



(\\X'-XY^^ + \\A'-AY^2,, 



111- 



I \Xl-Xl\\\Xl-XY^\A\-Al\^)ds 
Jo J 

< eeVf(\\x^ - x^%2 + \\v^ - y^||^2 + \\x^ - xY^^) , (4.10) 

for some constant 9e — 9q{C). Therefore, we can plug the estimates (4.5)-(4.10) 
back into (4.4), obtaining the inequality 

\\X^-XY^. < 9rVT{\\X'-XY^, + \\V'-V^%^ + W'XY^2), (4.11) 

for some constant 6*7 = 9j{C,K). On the other hand, thanks to the explicit 
nature of the functional differential equation (4.3), we can estimate ||V"^ — y^||j;«2 



15 



as follows: 



< E 



f{s,xly{x\v')s,c\M{x\v'))s) 



< TE 



- f{s, x^ y{x\ v^)s,c\M{x\ 
r f{s,xly{x\v%,c\M{x\v^))s) 

Jo 

- f{s, xly{x\ v^)s, cHm{x^ v^))s) 



ds 



ds 



< 3C^T\^E 

f-T 

e' 



IXl-X^l'ds 


+ E 


in 






.Jo 



[ \y{x\v'),-y{x^,v\ 

Jo 



ds 



< 3C^t(^\\x^ - x^f^2 + \\y{x\v^) - y{x^, v 

+ \\CHM{X\V^)) - C\M{X\ V^ml.), 



2m|2 



and by using the same arguments as for the estimates (4.5)-(4.10) this yields 
that, for a constant = 6%{C, K), 



,Vr(||xi- 



-x^%^ + \\v'-v 



^"|.2 + \\X^ - X^\\'^2 



)• (4.12) 



Hence, we can sum the inequalities (4.11) and (4.12), obtaining a constant 
6*9 = 9q{C, K) depending only on C and K such that 



< e,^/T{\\{X\V') - (X^ t/2)||2^^^^^ + \\{X\V') - (X^y2)||^. ,^.), 
which impUes that, for T > such that 9q\/T < 1/2, 



ll(x^vl)-(x^y2)ll^.,^. < 



- 1 



\\{x\v')-{x\v^)r^.^^.. 



<1 



Therefore, L is a contraction if T < ^{C, K) := A 1, and thus admits a 

unique fixed point {X, V). □ 

Remark 4.6. The extension of Theorem 4.5 to the case where £^ satisfies (LI) 
with respect to 0'^{[Q,T],W) = Jf'^{[0,T],W^) is problematic, as shown by 
the following counterexample borrowed from the theory of FBSDEs. Assume we 
have an augmented Brownian filtration, and consider the functional differential 
system 

'dXt = C''{M'^{X,V))tdWt, 
dVt = 0, 
Xo = Vo = 0, 
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where L'^ is the operator given by Ito's representation, and c6(x) := x + Wt- We 
thus obtain that V = 0, and the first equation can be rewritten as 



Assume it has an adapted solution X. We would then have Xt = E[XT + WT\J^t\ 
for all t > 0, which would lead in particular, for t = T, to the contradiction 
Wt = 0. 

Remark 4.7. While already quite general. Theorem 4.5 can be further ex- 
tended in a number of different ways. We present some possible extensions 
which only require slight modifications of the proof presented above: the details 

arc left to the reader. 

(i) First of all. Theorem 4.5 can be extended to any initial time r > 0: more 
exactly, consider functional differential systems on [r, T] of the form 

'dXt = nit, Xt,y{X, V)t, C\M{X, V))t)dt + a{t, Xt,y{X, V)t, C^MiX, V))t)dWt, 
< dVt = fit, Xt,yiX, V)t, C^MiX, V))t)dt, 

where r;,C G L'^{Tr). Then it is possible to prove that, under the same 
conditions of Theorem 4.5 and provided that T — t < i, there is a unique 
solution in S^'^i[T,T],W) x S^^i[T,T],W^). This remark will be essential 
in the next section, where we extend the solvability to any time interval. 

(m) It is not difficult to see that Theorem 4.5 also holds for operators 
and £^ of the form = • • • , Cl..,Cl.._^_-^, ...,£;.), i — 1,3, where £j 

satisfies Assumption (LI) with respect to {[0,T],W'i) for 1 < j < 

fcj, and with respect to Jf^{{0,T],W'^) for h + 1 < j < If. indeed, it 

suffices to consider the space ff'^{[0,T],Rp') := Yl''.l^S^'^{[0,T],Wi) x 

UlLki+i '^^([0,T],IRp5) and take the norm given by 



(m) The result of Theorem 4.5 can also be extended to other types of terminal 
condition. Indeed, assume that D([0,T],M"') denotes the space of all Re- 
valued cadlag functions on [0, T], and let $ : nxD([0, T], M") ^ satisfy 
the L°°-Lipschitz condition 



dXt = dM'^iX,0)t, Xo = 0. 




$(x) - $(x')| < C sup |xt - x;| P-a.s. Vx,x' e©([0,T],M") 



te[o,T] 



where C > 0. Then, the operators 



MiX,V)t := E[^iX) + VrlTt], y{X,V)t :=M{X,V)t - Vt 



satisfy estimates similar to those of Lemma 4.2. Therefore, the reader can 
easily verify that Theorem 4.5 remains valid if we substitute the operators 
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Jv[ and y in the system (3.3) byA^ and y. Moreover, the same conclusion 
remains true if $ satisfies the L^-Lipschitz condition 

|$(x) - $(x')| < C / |xt-x;|rfi P-a.s. Vx,x' e©([0,T],M"), 
Jo 

instead of the above L°°-Lipschitz condition. Two typical examples are 
the functionals $^(x) = supjgjg.T] 1^*1 ^'^'^ $^(x) = Xtdt, which are 
related to lookback and Asian options. 

{iv) Finally, it is possible to generalize the system (3.3) by substituting the 
equation for the component V by 

dVt = f{t,Xt,y{X,V)t,CHM{X,V))t)av{dt), Vo = 0, 

where, for some T > 0, ay is a non-random, finite Borel measure on [0, T] 
such that Q!y({0}) = 0. This is done by applying a change of integration 
order similar to the one discussed in (4.1). 

We conclude this section by observing that the flexibility in the choice of the 
operators opens the door to probabilistic interpretations for many classes 
of integro-partial differential equations, similarly to the well known non-linear 
Feynman-Kac formula for BSDEs. We do not attempt to investigate this prob- 
lem into more detail here, leaving it for future research. 

5 Solvability on arbitrary intervals 

After proving in Theorem 4.5 the local solvability of the fully coupled system 
(3.3), the next, more problematic step consists in showing the existence and 
uniqueness of solutions on arbitrarily large time intervals. In the case of simple, 
non-coupled Lipschitz functional differential equations, Liang et al. [15] showed 
that the existence and uniqueness of solutions can be obtained by imposing 
additional conditions on the operator £, while leaving the assumptions on the 
driver / and the terminal condition ^ unchanged. 

However, the situation is quite different for coupled systems: as it is now 
well known in the theory of classical FBSDEs, an extension of Theorem 4.5 to 
arbitrary intervals is possible only if we impose additional assumptions on the 
coefficients /i, cr, /, cf) of the fully coupled system (3.3), since the solution could 
explode for large time horizons (without going into more detail, we refer the 
reader to classical counterexamples for Markovian FBSDEs and related PDEs 
which can be found for instance in [18]). Several techniques have been adopted 
for classical FBSDEs to overcome this difficulty, as we mentioned in the Intro- 
duction. However, we prefer to apply another approach which in our opinion is 
the most natural in the case of functional differential systems. 

First of all, we briefly discuss the intuition. Similarly to [15], the first step 
consists in dividing the interval [0, T] into a finite number of subintervals Ij := 
[Tj_i,Tj], = To < ••• < Tjv = T: then, we solve the system separately 
on any subinterval, starting from the last one and going backward. There is, 
however, an important additional difficulty with respect to [15], which clarifies 
why additional assumptions on the coefficients are needed. For such an approach 
to work, we need in fact that the length of the subintervals Ij, on which the 



18 



system has to be solvable, can be bounded by below by a constant independent 
of J. We have seen that such a length depends on the Lipschitz constants C and 
K of the system, and we can observe that the only potential complication can 
arise from the terminal condition on each Ij : indeed, while the other coefficients 
ji, a, f and remain the same on all intervals Ij, the terminal condition of 
the subsystem on Ij is given by = y{X^+^ ,V^+^)t^ , where {X^+^,V^+^) is 
the local solution on Ij+i- To obtain the desired lower bound for the length for 
all Ij, it is therefore sufficient that S-' = Oj{X^^), where Oj{ui,x) is for each j 
Lipschitz continuous in x with respect to some constant C independent of j. 

Due to the strong coupling between the; two functional differential equations, 
this uniform Lipschitz continuity can be obtained only by studying the interplay 
between the two components X and V. However, this is very difficult without 
a concrete expression for the operators and is especially true in the case 
of non-local operators (namely, the fact that £(M) could still depend on the 
whole path of (Mt)tg[o ni^Y easily cause the explosion of the solution for large 
time intervals). These observations lead us to the conclusion that the solvability 
on arbitrary intervals of the fully coupled system (3.3) has to be treated on a 
case-by-case basis, by developing tailor-made techniques for each choice of the 
operators 

As an example, we study in the following the case where the system (3.3) is 
associated to a classical Brownian FBSDE. Hence, we will assume for the rest 
of the article that the filtration is generated by an m-dimcnsional Brownian 
motion {Wt)te[o,T] on {Q,J^,P), and that = C? = C? = £, where £ is the 
operator given by Ito's representation theorem. For notational simplicity, we 
will write Z(X. Y) = C{M{X, V)). 

In the literature of classical FBSDEs, it is critical to distinguish between the 
cases where the coefficients ^, a, f and (p are purely deterministic (i.e. they do 
not depend on to) or not. This is due to the fact that, in the purely deterministic 
case, one can exploit the well known connection between FBSDEs and parabolic 
PDEs. In the framework of functional differential equations, such an approach 
has been adopted by Liang et al. [14] in a very special case, where they essentially 
rely on the Lipschitz continuity of the solution of the corresponding parabolic 
PDEs (see for example Delarue [4,5]). 

We are, however, more interested in the case where the coefficients /i, a, 
f and 4> are allowed to be random. In the literature, there are several articles 
dedicated to the random case: in particular, Zhang [25,26] proved the solvability 
of classical FBSDEs by deriving imiform Lipschitz estimates with respect to the 
initial condition of the component X. In the following, after reformulating in 
our setting Zhang's results on uniform Lipschitz continuity, we briefly illustrate 
how we can use these results together with Theorem 4.5 to construct a solution 
to our coupled functional differential system. For the rest of this section, we 
will assume that n = 1, i.e. the component X is 1-dimensional, and that a does 
not depend on Z2, i.e. a = (7{t,x,y). Moreover, we assume that the coeSicients 
of (3.3) satisfy the following condition: 

Assumption (A2): The functions ^ : x [0, T] x R x R'' x M''^'" R, 
a : n X [0,T] X R X R'' ->■ M™, / : 17 x [0, T] x R x M'^ x R''^'" -)> R'^ and 
(f) : fl xR ^ R'^ satisfy Assumption (A2) if they satisfy Assumption (Al) with 
(A1.2) replaced by (A1.2'), and: 
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(A2.1) is uniformly Lipschitz in x with constant C. 
(A2.2) There is a constant 7 > such that 

Aliy) < -7|A?(j/)|, yy e {y' e R'' | \y'\ = 1}, 
for all y gR'^ with \y\ = 1, where 

d 

i=l 

+ dx(y{dy_ii)^y + {dyfj,)^y, 
Kl{y) := \dM' - lid.l^fyf + 2y^ d,y.{dyaf y, 

and where we assumed that all corresponding derivatives exist. 

The assumption on (j) is simply needed in order to differentiate between the 
Lipschitz constant of </> and those of the other coefficients. Observe moreover 
that we require stronger regularity conditions than those of Delarue [4,5], since 
the coefficients of the system (3.3) are random (however, a is allowed to be 
degenerate). Under these conditions, by deriving some clever estimates for 
linear FBSDEs, Zhang obtained the following result: 

Lemma 5.1 (Zhang [26]). Assume that the coefficients fx, u, f and (p satisfy 

Assumption (A2) with J — , and let £ denote the constant in Theorem 4-5. 
Let T < i, and for G R, i = 1,2, let (X% V^) denote the solution of 

dxi = ^i{t,xi,y{x\v%,z{x\v%)dt + a{t,xi,y{x\v%)dWt, 
dv; = f{t,xi,y{x\v%,z{x\v%)dt, 
XI = x\ = 0. 

Then, there is a constant qc, depending only on C, such that 
\y{X\V^)o-y{X\V^)a\<C\x^-x% 

where C := ^(|C'|2 + l)eecr _ 1 > q. 

The proof can be found in [26]. As already mentioned, Lemma 5.1 in con- 
junction with Theorem 4.5 allows to construct a solution to the system (3.3) on 
arbitrarily large time intervals. More exactly: 

Theorem 5.2. Assume that the coefficients jjL, a, f and 4> satisfy Assumption 
(A2). Then, for any T > 0, the fully coupled system (3.3) has a unique solution 

{X,V) in ^2([o,T],]R") x 5^'^{[Q,T],R'^). 

Proof. We assume without loss of generality that 7 = ^ in Assumption (A2), 
by changing C or 7 if necessary. Let C denote the constant in Lemma 5.1, 
and let £ = £{C) be the constant in Theorem 4.5. Wo consider a partition 
(To, Tiv) of [0, T] such that = Tq < • • • < Tat = T and < Ti-Ti_i < I for 
i = 1, - ■■ ,N, and set /, := [Tj_i, T,]. Moreover, to emphasize the dependence of 
the operators y, A4 on the terminal time and condition, we introduce a slightly 
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different notation and write, for all t G [0,Tj], 

M^'{^,V)t = E[^ + VT,\J't], 

y^^{^,v)t = E[^ + VT,\J't]-Vt. 

Z'^' V)t is then defined via the Ito representation of {M'^' {£,, V)t)t£[o,Ti] ■ 

The first step consists in constructing appropriate terminal conditions for 
all subintervals Ij. This is accomplished via a backward procedure. We set 
^jv := (l>, Cn := C, and for all a; e M, we consider the following system on In: 

■dX^-^ = Hit, Xf ■^ 3)^^ {0n(X^P, {On{X^P, V''nt)dt 

" dv,''''' ^ f{t,x^-'',y^-{9N{x^p,v^nt,z^''{0N{x^p,v^nt)dt, 

Xrj, = X, Vrj, ~ U. 

Since On has Lipschitz constant Cn < C, the system has for all a; S R a unique 
solution (X^'^, 1/^'^) by Theorem 4.5 and Remark 4.7 {i). We can thus define 

^"^"^ On-i(x) i^^"(0^(X^,^),F^'^)T„_,. 

It is then easy to check that 9n-i{x) is J^Tn-i -measurable for all a; € M. More- 
over, by Lemma 5.1, 9n-i is uniformly Lipschitz in x with constant 

Cn^i := ^(|Civ|2 + l)eec(T«-T„_i)_i. 

Since Cn-1 < C, we can iterate the same argument: for i = A'^ — 1, • • • , 2, 
we consider for all a: G M the solution (X*'^, V"''^) on li of the system: 

'dx*'" = i^{t,xi''',y'^'{ei{x'j^^),v''nt,z^'{eiix'^^),v''^)t)dt 

' dvr = fit,xi'''j^*iei{x^^^)lv'nt,z^'mx^^^),v'nt)^^^ 

which exists by Theorem 4.5 and Remark 4.7 {i). We then define 9i-i by 
^,_i(a;) :=y^^{ei{X'^^),V'nT,_,. 

6i-i{x) is thus -measurable for all x G R, and by Lemma 5.1, 6i-i is 

uniformly Lipschitz in x with constant Cj-i := (IC'iP + l)es<'(^*~'^'-i) — 1. 
Moreover, we can easily verify by induction that 

= ^{\CN\^ + l)eeoiTM-T,.,) < c. 

Now that wc have derived appropriate terminal conditions di, we can con- 
struct the solution on the whole interval [0, T] by a forward procedure. We set 
:= X, := 0. Then, for i = 1, • • • , N, we denote by {X\ the solution 
on /j of the system 

'dxi = XI, j)^- (e,(x^j, v')t, z'^- (eiix^T,), v%)dt 
' dv; = f{t,xi,y^^ie,ixirJ,v%,z'^*{0i{xi.j,v%)dt, 

Xi \ri — 1 T/i Tri — 1 
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which exists due to Theorem 4.5 and Remark 4.7 (i). We then set, for t G li, 

Xt:=Xi, Vt:=Vi. 
To prove that {X, V) solves (3.3) on [0,T], it suffices to check that, for t G li, 

y^^{eiix^j.j,v% = y{<f>{XTj,v)t, z'^^{6i{xir,),v% = z{4>{XT^),v)t. 

However, for z = 1, • • • ,N — 1 and t G [0, Tj], we have that 

= E[ei+i{XT,^,) + Vt,+, - Vt, + Vt. I - Vt 

by the construction of 9i. This gives by induction that y'^^ {Oi{Xlj,,),V^)t = 
y{eN{XT„).. V)t = y((t>{XTr,\V)t on /i for i = l, • • • , N. On the other hand, 
this impHcs that M'^^(e,{X\,^),V')t = M{(j){XT„),V)t on /, for all i. In 
other words, {M'-^ {9^[X}p.),V^)t)teii is the restriction to /j of the martingale 
(.M((^(Xt„), l^)t)te[o,T] and, due to the locality of the operator Z, this gives us 
that V% = Z{c^{Xt,-), V)t on 

This shows that {X,V) is a solution of (3.3) on [0,T], and the proof is 
concluded by observing that the uniqueness is a consequence of the uniqueness 
{X\V') on li. □ 

Wo conclude this article by briefly mentioning an extension of Zhang's re- 
sults [25,26] recently derived by Ma et al. [17]. Motivated by the connection 
between FBSDEs and PDEs in the deterministic case, the authors suggest that, 
for random coefficients, the solution can bo extended to arbitrary intervals by 
relying on the existence of a random field 9 (called decoupling field) such that 
y{X,V)i = 0{t,Xt). It is then sufficient to show that 6 is uniformly Lipschitz 
continuous. This can be obtained via the introduction of a backward stochas- 
tic Riccati equation of quadratic growth, under much weaker assumptions than 
Assumption (A2) (however, all processes have to be one-dimensional). 
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